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The lifetime of the B~ meson is measured using 272 exclusive B~ — > J/ijj(—} fj, + /j,~) n~ decays 
reconstructed in data from proton-antiproton collisions corresponding to an integrated luminosity of 
6.7 fb _1 recorded by the CDF II detector at the Fermilab Tevatron. The lifetime of the meson 
is measured to be t(B~) = 0.452 ± 0.048(stat) ± 0.027(syst) ps. This is the first measurement of 
the B^ meson lifetime in a fully-reconstructed hadronic channel, and it agrees with previous results 
and has comparable precision. 
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In the standard model the B~ meson is the only meson 
composed of two distinct heavy quarks. The B~ meson 
decay can be governed by the decay of the b or c spectator 
quarks or can proceed through the annihilation of the b 
and c quarks. Various theoretical techniques have been 
used to predict the B~ meson lifetime. An operator- 
product-expansion calculation [l| predicts a B~ meson 
lifetime in the range of 0.4 to 0.7 ps. A QCD sum rule 
approach 0] predicts the lifetime to be 0.48 ± 0.05 ps. 
Another approach Q, estimating the B~ meson lifetime 
by global fitting of the phenomenological parameters of 
all other heavy mesons, gives a result of 0.36 or 0.47 ps, 
depending on different choices of effective heavy-quark 
masses. 

The B~ meson lifetime was measured previously in 
semileptonic decays by CDF Q and DO Q. These 
measurements have an undetected neutrino in the final 
state and rely on the modeling of the B~ meson mo- 
mentum to account for unmeasured momentum of the 
neutrino. Therefore, a measurement of the B~ meson 
lifetime in a fully-reconstructed decay mode is desired 
since it does not suffer from this limitation. CDF is 
the first experiment to observe the fully-reconstructed 
B~ — > J/ip(-^_U + H~) 7T~ decay mode [6j and to measure 
the B~ mass [7[ . In this paper we present a lifetime mea- 
surement of the B~ meson using this decay mode. This 
measurement is made using data from pp collisions at a 
center-of-mass energy of 1.96 TeV recorded by the CDF 
II detector. The results are based on a data sample with 
an integrated luminosity of 6.7 fb _1 . 



The components of the CDF II detector [8J] most rel- 
evant for this analysis are the charged-particle tracking 
and muon identification systems. The tracking system 
is immersed in a uniform 1.4 T solenoidal magnetic field 
coaxial with the beam line, and consists of single- and 
double-sided silicon detectors @ and a 96-layer open-cell 
drift chamber (COT) The muon system is used 

to identify the J lip l-i + M~ decay. Two sets of muon 
chambers [lj, [12] are used to cover different pseudora- 
pidity regions. 

A three-level event-selection system (trigger) is used 
to collect events enriched in Jjip — > fi + \T~ decays [|[. 
The event reconstruction starts by combining two muon 
candidates to form a J/ip candidate. The trigger require- 
ments are confirmed by selecting events that contain two 
oppositely charged muon candidates, each with matching 
COT and muon chamber tracks. The muon-pair mass is 
required to be within 80 MeV/c 2 of the world-average 
J/tp mass [HI], where the mass resolution is 13 MeV/c 2 . 

Both J/ip K~ and J/ip tt~ combinations are recon- 
structed in this analysis. The large B~ —> J/tp K~ sam- 
ple is used as a reference decay for B~ — )• J /if) n~ . These 
final states are identified by assigning the K ~ or ir~ mass 
to other reconstructed tracks not used in the J/ip can- 
didates and forming B~ or B~ candidates. Each three- 
track combination must satisfy a kinematic fit in which 
the three tracks are required to originate from a com- 
mon decay point, and the invariant mass of the muon 
pair is constrained to the world-average J/ip mass [13l |. 
The K~ and tt~ candidate track is named the third track 



h~ . A minimal selection is made on kinematic quantities 
after the constrained fit including p T {hr) > 1.7 GeV/c 
and pr{B) > 5 GeV/c, where pr is the momentum com- 
ponent transverse to the beam line, and B refers to a 
J ftp h~ candidate. The selection criteria for the B can- 
didates are listed in Table U and discussed below. 

TABLE I. Selection requirements. 



Selection variable 



Requirement 



pr(h ) 
Pt(B) 

p(x 2 ) 

\d{B)\/<J d(B) 

Pt 

Ib 
o m {B) 

ct 
o-ct(B) 



> 2.0 GeV/c 

> 6.5 GeV/c 

> 0.1% 
< 2.0 

< 0.2 radians 

> 0.6 

< 40 MeV/c 2 

> 80 ^m 

< max[35, 65 - 3 p T (S)(GeV/c)] fim 




5.2 5.3 

Mass(J/>|/ K") [GeV/c 2 ] 



FIG. 1. Invariant-mass distribution of J/tj) K~ candidates. 
The hatched areas are the sideband regions and the signal 
region lies between them. 



The h~ and B candidates are required to have a min- 
imum pt to suppress combinatorial background events. 
We reject events with poorly defined decay points by im- 
posing a lower threshold to the chi-square probability 
P(x 2 ) °f t ne constrained fit used to reconstruct the B 
candidates. We select B candidates that originate from 
the primary interaction point by requiring a small impact 
parameter d(B) (transverse distance of closest approach 
to the beam line) in units of its uncertainty pmb) and a 
small angle 0t between Lt and pr(B), where Lt is the 
transverse displacement vector from the primary inter- 
action point to the .B-decay point, and pt{B) denotes a 
vector in the transverse plane along the momentum di- 
rection of the B candidate. The isolation Ib of the B 
candidate is defined as Is = p(B)/{p{B) + |X)iPil)> 
where ^2 i Pi is the sum of momenta over all other re- 
constructed tracks not used in the J/ip h~ combination 
within \/(Arj) 2 + (A0) 2 < 0.7, and A-q and A0 are the 
differences in pseudorapidities and azimuthal angles of 
tracks relative to p{B). We also suppress the promptly 
produced combinatorial background by rejecting candi- 
dates with small ct, where ct is the decay length of the 
B candidate determined by 



ct 



L T -pt(B] 



c m(B) 
\pt(B)\ 2 



(1) 



and m(B) is the reconstructed mass of the B candidate. 
Requirements on a m (B) and a c t(B) are made to reject 
poorly reconstructed events, where a m (B) and a c t{B) 
are the associated uncertainties from the kinematic fit 
of m(B) and ct(B), respectively. The optimization of 
the selection requirements is obtained by maximizing the 
quantity S/y/(S + B) where the background B is esti- 
mated from the mass sidebands in data and the signal S 
is estimated from the signal-region data after subtracting 
background contributions as calculated from sidebands. 



Because the a c t(B) distribution depends otipt(B), we 
vary the requirement on a c t(B) as a function of pt(B). 
For candidates with pr{B) < 10 GeV/c, we require 
o'ct(B) < (65 — 3 pr{B)) jttm for pr{B) measured in 
GeV/c, and a ct {B) < 35 nm for p T {B) > 10 GeV/c. 
This py-dependent requirement on a ct (B) is chosen to be 
highly efficient for preserving signal while reducing com- 
binatorial background and leads to no measurable biases. 

The resulting B~ mass distribution is shown in Fig. [T] 
The signal region lies between two background sideband 
regions and has 46 280 B~ candidates. The two side- 
band regions consist of a lower sideband from 5.18 to 
5.23 GeV/c 2 and an upper sideband from 5.33 to 5.38 
GeV/c 2 , as shown in the hatched areas. 

The B~ — > J/ip tt~ candidates are formed from the 
same parent sample as the B~ r- J/tjj K~ candidates 
where the only change to the reconstruction is to assign 
the pion mass to the third track. We then select events 
for further analysis using the selections in Table [J] The 
reconstructed mass distribution for the B~ candidates is 
shown in Fig. [5] The signal region lies between two back- 
ground sideband regions and has 1496 B~ candidates. 
The two sideband regions of B~ candidates consist of 
a lower sideband from 6.16 to 6.21 GeV/c 2 and an up- 
per sideband from 6.33 to 6.60 GeV/c 2 , as shown in the 
hatched areas. The lower sideband is narrow to avoid 
contamination from semileptonic B~ decays where the 
lepton is misidcntified as a pion. 

We generate Monte Carlo simulations for 
B~ — > J/ip K~ and B~ —> J/ip ir~ decays to study 
the efficiency of the selection criteria as a function of 
decay length. The Fixed-Order plus Next-to-Leading 
Logarithms (FONLL) pr spectrum [l4| is used for 
the B~ production spectrum. We use the calculation 
of Ref [15| as the spectrum for B~ simulation. In 
comparing the FONLL px spectrum for B~ production 
with the experimental data, reasonable consistency is 
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simulation are also pT-dependent, and the px thresh- 
old remains the same; the only change is to require 
a ct (B) < (45-2 p T {B)) ^m for p T (B) < 10 GeV/c, 




FIG. 2. Invariant-mass distribution of J/ip n~ candidates. 
The hatched areas are the sideband regions and the signal re- 
gion lies between them. The fit result is overlaid in the signal 
region, as well as the signal and background components. 



i 1 1 1 1 1 1 1 1 1 1 1 1 

0.02 0.04 0.06 0.08 0.1 0.12 

ct(JA|/ if) [cm] 



FIG. 4. Comparison of efficiency for B~ — ► J/ip K~ obtained 
from data and the fit result from simulation. Also shown is 
the fit result for B~ — > J/ip %~ simulation. 



Data 
-— Simulation 



0.004 
c a (J/y K ) [cm] 



FIG. 3. The a c t distribution of J/ip K candidates obtained 
from the simulation is compared with data. 



observed for px > 6 GeV/c, and any residual discrepancy 
gives a negligible systematic uncertainty. To further 
validate the B~ — >• J/ip K~ simulation, we compare the 
distributions of the selection variables listed in Table |l] 
for experimental data and simulation. Generally, good 
agreement is observed for all selection variables except 
for a c t(B), whose comparison is shown in Fig. [3J The 
disagreement in the a c t(B) distribution arises from 
mis-modeling of the silicon tracking detector in the 
simulation, giving a smaller a ct (B) compared with 
experimental data for a given px(B). Consequently, the 
selection requirement made on a ct (B) for simulation is 
tuned in order to allow the same efficiency as in the 
experimental data. These <r ct {B) selection values for 



and a ct {B) < 25 pm for p T {B) > 10 GeV/c. The 
systematic uncertainty associated with the tuning will 
be discussed later. 

The efficiency of the selection criteria is found by com- 
paring the decay-length distribution after applying the 
selection in Table U to that obtained from the mini- 
mal selection which requires px(h~) > 1.7 GeV/c and 
pr(B) > 5 GeV/c. The efficiency determined from sim- 
ulation is then fit to a function of the form 



e(ci) = C 



1 — exp 



ct 



(2) 



where C, a, b are parameters to be fit. Figure H] shows 
the efficiency determined from B~ — > J/ip K~ experi- 
mental data as well as the fit result from simulation. 
The parameter C in Eq. ([2]) is not necessary in the life- 
time fit because only the relative shape of the efficiency 
function matters. The requirement on /3t leads to an 
efficiency that is not constant as a function of decay 
length. This variable is very effective in rejecting back- 
ground events, especially for events with small ct. The 
good agreement between the simulated efficiency and the 
data-determined efficiency supports the use of this ap- 
proach in the ct-dependent efficiency. The efficiency for 
the B~ — ► J/ip 7T~ decay as a function of decay time de- 
termined from simulation is fit and also shown in Fig. 2) 
We use a maximum log-likelihood simultaneous fit to 
the unbinned mass and decay-length distributions of the 
B~ candidates. The likelihood function consists of signal 
and background parts, and each part has a mass term and 
a decay length term. The log-likelihood function is 
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ln£ = ln ifa M s{m) T s {cti) + (1 - /.) M b (m t ) T b {cti)\ , (3) 

i 

i 



where f s is the signal fraction, m, and cti are the re- 
constructed mass and decay length for event i. M s (rrii) 
and T s (cti) are the normalized probability density func- 
tions for mass and decay length of the signal model, and 
Mb(rn,i) and T^icti) are the corresponding functions of 
the background model. The signal mass model M s (mi) 
is described by a Gaussian distribution with mean mo 
and width o m , whose values are determined by the fit. 
The signal decay length model T s (cti) is an exponential 
distribution with characteristic lifetime t, smeared by the 
detector resolution and multiplied by the efficiency func- 
tion given in Eq. @. The detector resolution, which is 
modeled as a Gaussian distribution centered at zero with 
a width of 20 /an, is chosen to be consistent with cali- 
bration using promptly decaying background events [l6[ . 
The background mass model Mb(nii) is described by a 
linear distribution, and If, (cti) is described by a linear 
combination of three exponential distributions. 

A two-step process is used to extract the lifetime of 
the B~ meson. The first step includes the efficiency 
fit and the sideband fit. The efficiency fit is performed 
on the simulated events using Eq. @, and the result is 
shown in Fig. |4j The sideband fit consists of two sepa- 
rate fits in the sideband regions. The first fit determines 
the background mass model parameters, and the second 
fit determines the background decay-time model parame- 
ters. The signal region is fit with the efficiency and back- 
ground parameters Gaussian-constrained to the result of 
the earlier fits. The signal fraction and the signal model 
parameters are allowed to float freely in maximizing the 
log-likelihood function given in Eq. ([3|). 

To validate this fitting technique, the fit is first applied 
to the B~ candidates shown in Fig. [1] to extract the B~ 
meson lifetime. The B~ lifetime which is found to be 
t(B~) = 1.637 ± O.OfO (stat) ps is in agreement with 
the world-average value of 1 .641 ± 0.008 ps pj]. The B~ 
signal yield returned from the fit is 43 308 ± 171. The 
fit is then applied to the B~ candidates shown in Fig. [2j 
resulting in a lifetime of r(B~) = 0.452 ± 0.048 (stat) ps 
for the B~ meson, which is taken as our central result, 
with a B~ signal yield of 272 ± 61 (stat) candidates. The 
B~ meson mass of 6274.6 ± 2.9 (stat) MeV/c 2 returned 
from the fit is in good agreement with the previous CDF 
determination Q ■ Figures [2] and [5] show the distribution 
of the J/tp ir~ mass and decay length. The fit projections 
are overlaid. 

We have considered several sources of systematic un- 
certainty and evaluated their contributions. To evaluate 
possible systematic uncertainties with the models in the 
likelihood function, we generate 400 simulated data sam- 
ples (pseudoexperiments) whose distributions are based 
on the fit results determined by the experimental data. 
These pseudoexperiments are then fit with the default 




ct(JA|/ jc") [cm] 



FIG. 5. Decay-length distribution of J/ip ir~ candidates. The 
fit projection, along individual contributions from signal and 
background, is overlaid. 



and alternate models separately. The distributions of the 
sample-by-sample lifetime differences between different 
models are obtained and compared with the differences 
observed in experimental data. To assess the effect of the 
choice of the linear model for the mass-fit background, 
we compare to the result of a fit using a bilinear model 
that allows the background distribution to have different 
slopes at masses lower and higher than the B~ pole mass, 
with the constraint that these two distributions intersect 
at the fit B~ mass value. The fit lifetime with this bi- 
linear model has a shift of -0.009 ps compared with the 
default linear model. The pseudoexperiments suggests 
up to a 0.017 ps difference from this variation. We con- 
clude that the shift between the data fits is consistent 
with the spread among the pseudoexperiments, and we 
use that larger difference as the systematic uncertainty 
from the background mass model. 

To assess a possible systematic uncertainty due to 
the modeling of the long tail in the background decay- 
length distribution, we test an alternate model of the 
background decay time which uses a linear distribu- 
tion to replace the component with the largest charac- 
teristic lifetime in the three exponential distributions. 
This variation gives a lifetime result that changes by 
-0.0007 ps compared with the default background decay- 
time model. However, fit results from pseudoexperiments 
suggest the difference between these two models could be 
0.013 ps, which is included as the systematic uncertainty 
due to the choice of the background decay-time model. 

The signal decay-time model includes the efficiency de- 
termined from the simulation. We have performed sev- 
eral studies to estimate the associated systematic un- 
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certainty. First, the fit is repeated using an efficiency 
function obtained without tuning the a ct (B) difference 
between data and simulation. The difference in the es- 
timated lifetime is 0.003 ps. Second, the efficiency func- 
tion is shifted toward lower and higher decay length by 
20 /im to account for a possible uncertainty in determin- 
ing the efficiency function parameters; this 20 fim shift 
is equivalent to three standard deviations of the param- 
eter a in Eq. ([2]). This variation gives a difference of 
-0.01 (+0.007) ps for shifting toward lower (higher) de- 
cay lengths. The distribution of the difference between 
the resulting lifetimes in the pseudoexperiments is fit by 
a Gaussian distribution that centers at -0.006 (0.004) ps 
with a width of 0.002 (0.001) ps for shifting toward lower 
(higher) decay lengths. Third, the systematic uncer- 
tainty associated with the B~ production spectrum has 
been assessed. We vary the relative fraction of different 
contributions to the production spectrum; the difference 
in the corresponding efficiency is negligible and no sys- 
tematic uncertainty is assigned to it. Finally, to further 
study the systematic uncertainty associated with the pro- 
duction spectrum, we use the efficiency parameters ob- 
tained from the B~ — > J ftp K~ simulation. Since the B~ 
production spectrum is quite different from that of B~ , 
the fit lifetime difference of 0.007 ps indicates that the 
production spectrum does not contribute significantly to 
the systematic uncertainty. Thus, the total systematic 
uncertainty associated with the signal decay-time model 
is taken to be 0.010 ps. 

Correlations between the lifetime and other parame- 
ters of the analysis are considered as possible systematic 
uncertainties. The list of parameters includes the min- 
imum and maximum decay length for events in the fi- 
nal fit, adding a parameter to the efficiency model, small 
variations in the sideband definitions, small modifications 
in the selection requirements, the use of an alternate fit 
procedure which fits the sideband and the signal regions 
simultaneously, the mass resolution in the signal model, 
the background fraction, and the three terms describing 
the exponentials in the background decay time model. 
No systematic effect was found to significantly exceed 
the variations expected from statistical uncertainties. We 
assign an additional uncertainty of 0.010 ps as a conser- 
vative approach to account for possible small systematic 
effects. 

The systematic uncertainty due to tracking detector 
misalignments is evaluated by generating simulated sam- 
ples with radial displacements of individual sensors as 
well as translation and rotation of the silicon detec- 
tor relative to the COT pjj. A systematic uncertainty 
of 0.007 ps is assigned to the misalignment based on 
these simulated samples. The systematic uncertainty 
from the signal mass model is evaluated by including 
a contribution to the total B~ signal yield from the 
Cabibbo-suppressed decay B~ — > J/ijj K~ in the signal 
mass shape. The Cabibbo-suppressed contribution is 
fixed to be 5% of the total B~ signal yield as determined 
from the Cabibbo angle. This effect results in a 0.003 ps 



variation. 

The systematic uncertainty from the fitting technique 
itself is tested by generating pseudoexperiments, and 
comparing the fit lifetimes with the input lifetime. The 
bias on the lifetime returned by the fit is found to be 
no greater than 0.003 ps which we take as systematic 
uncertainty. Table |TT] summarizes the systematic uncer- 
tainties, which arc added in quadrature to determine the 
total systematic uncertainty. 



TABLE II. Summary of systematic uncertainties. 



Source 


Uncertainty [ps] 


Background mass model 


0.017 


Background decay-time model 


0.013 


Signal decay-time model 


0.010 


Correlation 


0.010 


Misalignment 


0.007 


Signal mass model 


0.003 


Fitting technique 


0.003 


Total 


0.027 



Given that the efficiency is not uniform over decay 
length, our result relies on the accuracy of the simula- 
tion in determining the efficiency. We check our result 
by measuring the B~ lifetime using a different set of se- 
lection criteria, each of which has uniform efficiency in 
the B~ — > J/ip K~ decay. The most important differ- 
ences between these selection criteria and those listed in 
Table U are removing the /3t requirement and using a 
larger minimum ct requirement. The alternate selection 
criteria gives 6538 B~ candidates between 6.0 and 6.6 
GeV/c 2 which is roughly the same number as obtained 
from the selections in Table Q] (6368 candidates), while 
only 2578 candidates are common to both samples. 

The consistency check also uses an unbinncd maxi- 
mum log-likelihood fit to extract the B~ meson lifetime. 
The signal mass model consists of a Gaussian distribu- 
tion centered at the B~ meson mass and a Cabibbo- 
suppressed B~ — > J/tfi K~ Gaussian distribution cen- 
tered at 60 MeV/c 2 below the B~ meson mass with a 30 
MeV/c 2 width. The signal decay-time model is an expo- 
nential distribution convoluted with the detector resolu- 
tion. The background mass model is described by a linear 
distribution, and the background decay-time model con- 
sists of two prompt Gaussian distributions, two positive 
exponential distributions, and one negative exponential 
distribution. 

A similar two-step fit is used in the consistency 
check. The first step is to determine the background 
parameters from the sideband fit, where the sideband 
is defined as the J '/ip n~ invariant-mass region between 
6.4 and 6.5 GeV/c 2 . The sideband fit is performed on 
events with decay length between -1000 /im and 1000 
fj,m and the resulting background parameters are fixed 
in the second step. In the second step we fit events in 
the signal region between 6.16 and 6.36 GeV/c 2 , and 
only the signal fraction and signal model parameters 
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are allowed to float. The consistency check is first 
performed on the B~ candidates. The fit result finds the 
B~ lifetime to be r = 1.647 ± 0.020 (stat) ps, 
which agrees with the world-average value of 
1.641 ± 0.008 ps QJ]. The consistency check is 
then applied to the B~ candidates, giving a B~ meson 
lifetime of r = 0.450 ± 0.053 (stat) ps which is consis- 
tent with our central value of 0.452 ± 0.048 (stat) ps. 
The B~ signal yield from the consistency check is 
308 ± 39 (stat) which is compared with 272 ± 61 (stat) 
from the central result. The total systematic uncertainty 
in the consistency check is 0.033 ps where the largest 
uncertainty of 0.027 ps comes from the background 
decay-time model. Thus, we conclude that our central 
result obtained from the ci-dependent efficiency is 
reliable. 

In conclusion, we have made the first measurement of 
the B~ meson lifetime in a fully-reconstructed hadronic 
decay mode. Using the B~ — > J/ip tt~ decay channel, 
the lifetime of the B~ meson is measured to be t = 
0.452 ± 0.048(stat) ± 0.027(syst) ps. This result is 
consistent with the most recent result from the DO col- 



laboration using semileptonic decay channels, r = 

0.448 iao3e;( stat ) ± °- 032 ( s y st ) P s > and has comparable 
precision. The result also agrees with theoretical cal- 
culations in which the decay width of the B~ meson is 
dominated by the decay of the charm quark. 
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